甲烷部分氧化过程中强制振荡的动力学Monte Carlo模拟

为了有效控制动力学振荡行为以提高反应的转化率, 利用Monte Carlo 方法研究了甲烷催化部分氧化过程中的强制振荡行为, 探讨了原料气周期性变化对振荡动力学和转化率的影响. 研究表明原料气周期性变化不仅可以有效调控振荡的动力学行为, 产生如短周期振荡和双峰振荡等特殊动力学过程, 而且还可以提高反应的转化率. 当强制过程的周期从T/3 增大到2T(T为自发振荡过程的平均周期), 振荡过程从短周期振荡变化为双峰振荡. 对反应过程中CO的转化率进行了计算, 结果表明原料气周期性变化可以有效提高反应的转化率.振荡动力学的改变和转化率的提高主要是因为强制振荡过程使得催化剂表面发生了从氧化态向还原态的转变.     

不同孪晶界密度银纳米线拉伸形变行为的分子动力学模拟

采用分子动力学方法模拟了不同孪晶界密度银纳米线的拉伸形变行为, 分析了孪晶界密度对多晶银纳米线屈服强度、弹性模量和塑性变形机理的影响. 在弹性形变区域, 孪晶界的存在对杨氏模量变化的作用不明显. 在塑性形变阶段, 首先从表面边缘开始产生位错成核, 然后延伸并受阻于孪晶界. 在进一步拉伸载荷作用下, 孪晶界将作为位错源产生新的位错. 模拟结果表明, 银纳米线的强度与孪晶界和晶粒的尺寸有关. 孪晶界密度较小(即晶粒的长径比大于1)时, 此纳米线的屈服应力比单晶纳米线还要小, 只有当孪晶界密度较大时(即晶粒的长径比小于1), 孪晶界使得纳米线得到强化. 综合分析了孪晶界和晶粒尺寸对银纳米线的影响, 为构建高强度金属纳米线打下基础. 最后讨论了温度和拉伸速度对孪晶纳米线屈服应力所产生的影响, 随着温度的升高, 孪晶纳米线与单晶纳米线的屈服应力差先升高后趋于稳定; 当拉伸速度逐渐增大, 孪晶纳米线与单晶纳米线的屈服应力差先稳定后增大.     

分子动力学模拟研究[Bmim][PF6]+水+醇三元体系的微观结构和分子间相互作用

研究离子液体体系的微观结构和分子间相互作用具有重要意义. 本文对1-丁基-3-甲基咪唑六氟磷酸盐([Bmim][PF₆])+水+乙醇和[Bmim][PF₆]+水+异丙醇三元体系进行了分子模拟研究, 计算了径向分布函数和不同组成的水-醇混合溶剂与离子液体阴阳离子间的相互作用能, 并将其分解为库仑相互作用能和Lennard-Jones(LJ)势能. 在此基础上, 研究了溶液体系的微观结构、分子间相互作用和相行为. 结果表明, 水倾向于与离子液体阴离子和阳离子极性部分作用, 醇倾向于与阴离子和阳离子非极性部分作用; 库仑力主导阴离子-溶剂相互作用, 色散力主导阳离子-溶剂相互作用, 阴阳离子的缔合状态对色散力影响较小, 对库仑力的影响非常显著.     

大环多胺金属模拟酶催化水解梭曼机理研究

为解决生物酶在战场应用环境中易失活且环境耐受性差的问题,我们合成了大环多胺双锌模拟酶配合物-2,6-双{[双(2-羟乙基)氨基]甲基}-4-甲基酚(Zn L2).Zn L2对GD(Soman,梭曼)具有较高的催化活性,能够实现多次催化循环,并且在Zn L2与GD化学计量比为(0.56∶1)的条件下水解率为100%.使用高效液相色谱-质谱与核磁共振确定催化产物为甲氟膦酸酯、甲氟膦酸、甲基膦酸频哪酯与甲基膦酸,说明Zn L2水解GD为双平行反应,反应途径分别为GD的P—F与P—O键的断裂.建立在实验的基础上,使用密度泛函理论(DFT)的B3LYP方法,从P—O与P—F键断裂这两个平行反应方向对催化过程进行计算.计算结果表明:Zn L2上Ot末端醇盐作为一个碱试剂去激活水分子,使其成为羟基并作为亲核试剂进攻中心磷原子,形成一个五配位的磷酰中间体,通过进一步亲核进攻水解中间体形成最终P—O与P—F键断裂产物.Zn L2催化GD水解的活化能为分别为5.6 kcal/mol(P—F键断裂)和11.5 kcal/mol(P—O键断裂),说明P—O键断裂过程是GD酶促水解反应的速控步骤.     

Simulation Methods for Functional Materials

Functional materials embrace a broad area, ranging from functional information materials to special polymers, from special chemicals for printing to matenais used in making paints. Inasmuch as most tunctional materials are particulate, the present contribution is considered pertinent to the present FORUM.     

随机结构正交展开分析的Ritz动力聚缩法

针对随机结果正交展开理论计算上的弱点 ,本文在分析扩阶矩阵特性的基础上 ,于 Ritz模态向量子空间中对扩阶方程实现动力聚缩 ,大大提高了正交展开理论对实际工程问题的分析能力。分析实例表明 :即使结构参数具有很大变异性 (如δ =0 .4 )时 ,该算法依然能理想地与 Monte Carlo法模拟结果相吻合 ,计算时间则远远小于 Monte Carlo模拟法。同时 ,分析例证再一次强化了在结构动力分析中考虑结构参数随机性的必要性     

平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体离散元法的增量平衡方程；而后推导了多刚体离散元的弹塑性弹簧系数矩阵，建立了多刚体离散元内力屈服面塑性铰法的增量求解格式，成功地进行了平面钢框架的弹塑性、大位移极限承载力分析。计算结果与其他数值方法或实验结果吻合良好，显示了多刚体离散元方法进行结构极限承载力分析这一复杂问题的优越性     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

有限元分析中硬化弹塑性本构方程积分算法

基于增量段上的应变关于塑性拉氏乘子的变化率为常矢量的假定,导出一种精确有效的、对有限元分析中线性混合硬化弹塑性本构方程的积分算法。以精确算例和等误差分布图的形式检验了本文方法以及其他两种常用方法的精度。